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Abstract We present the symplectic algorithm in the Lagrangian formalism for the Hamiltonian systems by virtue 
of the noncommutative differential calculus with respect to the discrete time and the Euler-Lagrange cohomological 
concepts. We also show that the trapezoidal integrator is symplectic in certain sense. 
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1 Introduction 

It is well known that the symplectic algorithml^^^l 
for the finite dimensional Hamiltonian systems are very 

. powerful and successful in numerical calculations in com- 
parison with other various non-symplectic computational 
schemes since the symplectic schemes preserve the sym- 

' plectic structure in certain sense. On the other hand, the 

. Lagrangian formalism is quite useful for the Hamiltonian 
systems. Since both are important at least in the equal 

' footing, it should not be useless to establish the symplec- 

■ tic algorithms in the Lagrangian formalism. As a matter 
of fact, the Lagrangian formalism is more or less earlier to 

I be generalized to the infinite dimensional systems such as 
classical field theory. 

In this note we present the symplectic geometry and 

, symplectic algorithm in the Lagrangian formalism in ad- 
dition to the Hamiltonian one for the finite dimensional 
Hamiltonian systems with the help of the Euler-Lagrange 
(EL) cohomological concepts introduced very recently by 

, the authors in Ref. [5]. 

I In the course of numerical calculation, the "time" 
' t (z R is always discretized, usually with equal spacing 
h = At, 

teR^te T = {(4-, ife+i =tk + h, keZ)}. 

It is well known that the differences of functions on T with 
respect to T do not obey the Leibniz law. In order to 
show that the symplectic structure at different moments 
tk is preserved, some well-established differential calculus 
should be employed. This implies that some noncommu- 
tative differential calculus (NCDC) on T and the function 
space on it should be used even for the well-established 
symplectic algorithms. In this note we employ this simple 
We also show that the trapezoidal integrator 
is symplectic in certain sense. Finally, we end with some 



2 The Necessary and Sufficient Condition for 
Symplectic Preserving Law 

In this section, we first recall some well-known con- 
tents in the Lagrangian formalism for the finite dimen- 
sional Hamiltonian systems. We employ the ordinary cal- 
culus to show that the symplectic structure is preserved by 
introducing the EL cohomological concepts such as the EL 
one-forms, the null EL one-form, the coboundary EL one- 
form and the EL condition and so forth. It is important 
to emphasize that the symplectic structure-preserving is 
in the function space in general rather than in the solu- 
tion space of the EL equation only. The reason will be 
explained later. 

Let time t e i? be the base manifold, M = M"- the 
configuration space on t, 

the (canonical) coordinates on it, T the transport, TM 
the tangent bundle of M with coordinates 

F{TM) the function space on TM. 

The Lagrangian of the systems under consideration is 
L{q^ ^ '?■') with the well-known EL equation from the vari- 
ational principle 



Let us introduce the EL one-form^^l 

It is clear that the EL equation is given by the null EL 
one-form, 

E{q\q')^Q, 

which is a special case of the coboundary EL one-forms 
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where a{q^,q^) is an arbitrary function of {q^,q^) in the 
function space F{TM). 

Taking the exterior derivative d of the Lagrangian, we 

get 



dL{q\q^)=E{q\q^) + —0, 
where 6 is the canonical one-form defined by 



(4) 



(5) 



Making use of nilpotency of d, 

d^L{q,q)=0, 

it follows that iff the EL one-form is closed with respect 
to d, i.e., 

dE{q\q^) = 0, (6) 

which is called the EL condition, the symplectic conser- 
vation law with respect to t holds, 

d 



dt 



LJ = 0, 



(7) 



where the symplectic structure u) is given by 

It is important to note that although the null EL one- 
form and the coboundary EL one-forms satisfy the EL 
condition and they cohomologically trivial, this does not 
mean that the closed EL one- forms are always exact. As 
a matter of fact, the equation (8) shows that the EL one- 
form is not always exact since the canonical one-form is 
not trivial in general. In addition, it is also important 
to note that the q'^{t), i = 1, . . . ,n in the EL condition 
are still in the function space in general rather than in 
the solution space of the equation only. This means that 
the symplectic two-form oj is conserved with respect to t 
with the closed EL condition in general rather than in the 
solution space only. 

In order to transfer to the Hamiltonian formalism, we 
introduce the canonical momentum 



- ^ 

~ dqi ' 



(9) 



and take a Legendre transformation to get the Hamilto- 
nian function 



H{q\pj)=pkq''-L{q\q^). 



(10) 



Then the EL equation becomes the canonical equations as 
follows: 



q' = 



dH 



p) = — 



dH 



(11) 



It is clear that a pair of the EL one-forms should be in- 
troduced now, 

E2{q\Pj) = -(^Pj+^^dq^. (12) 

In terms of z'^ = (q'% . . . ,Pn), the canonical 

equations and the EL one-form become 

z = J-^W^H, (13) 
noa2E{z) = dz^{Jz - V^if) . (14) 

Now it is straightforward to show that the symplectic 
structure-preserving law 

d 

di^ = '' 

u = dz"^ A Jdz (15) 
holds if and only if the (closed) EL condition is satisfied 
dE{z) = 0. (16) 

3 The Necessary and Sufficient Condition for 
Discrete Symplectic Preserving Law 

Now we consider the symplectic structure-preserving 
of symplectic algorithm in the Lagrangian formalism. As 
was mentioned above, in the course of numerical calcula- 
tion, the "time" i G i? is discretized with equal spacing 
h = At, 

t€R^teT = {{tk,tk+i=tk + h, keZ)}. (17) 
At the moment tk, the coordinates of the space 
G = {■ ■ ■ X ■ ■ ■ X ■ ■ ■} 
are denoted by q^'^^ , the symplectic structure by 

J'') = dq'i''^^ A dg^'^) , (18) 
and g^^^ now is the (forward-) difference of q^''\ 

A,g(") : = = qi''^ = ^(gC^+i) - g^^)) . (19) 

Now the EL equation becomes the difference discrete 
Euler-Lagrange (DEL) equation which can be derived 
from the difference discrete variational principle, 



(ft) 



dL 



D A 



dL 



(fe-1) 
D 

:(fc-i) 



= 0. 



(20) 
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Now we consider the difference discrete symplectic structure and its preserving property. Taking the exterior 
derivative d on T*(M^^), we get 

d4^) = ^d/'=) + |^|d,f). 
By means of the modified Leibniz law with respect to Af and introducing the DEL one-form, 

we have 

dL^') = E^^^ + M^^^ , (22) 



where 6^0^ is the discrete canonical one-form, 



oP-^d/'K (23) 
then there exists the following discrete symplectic two- form on T*{Mj.^), 

.(f) = ^e^^ = ,fff" , d,^-^''-^^ A d,^^'^^ + ^?#^d,^(. - 1) A d,^^'^^ . (24) 

(k) 

It is easy to see that the null DEL one- form Ejj' = gives rise to the DEL equation and it is a special case of the 
coboundary DEL one-forms 

E^^^ = dai'\q^''\qf'), (25) 

where an\q'^^''\ qt^ ^) is an arbitrary function of {q^^^\qi^ "*). 

Finally, due to the nilpotency of d on T*{Mj,^) it is easy to prove from Eq. (22) that iff the DEL one- form satisfies 
what is called the DEL condition 

AE^^^ = , (26) 

i.e., the DEL one-form is closed, the discrete (difference) symplectic structure-preserving law holds, 

Atw^f ) = . (27) 

Similar to the continuous case, the closed DEL one-forms are not always exact and this difference discrete symplectic 
structure-preserving law is held in function space in general rather than in solution space only. 
Let us consider the following DEL equation 

q{k) _ 2qik+i) ^ ^(fc+2) ^ _^2^( (fe+i)) _ (28) 

dq 

Introducing the DEL one-forml^l 

E^^+'^ : = d(g^('=+i)){g« - 2g('=+i) + q^'+^^ - h^^{q^'^+'^)} , (29) 

the null DEL one-form is corresponding to the DEL equation and the DEL condition directly gives rise to 

dqnk+^) A dg(^+i) = d^^f'^+i) A dg^ . 

It follows that 

A,a;W = ^[dql^'^'^ A dg^'^+i) - dqf^"^ A dq^'^^] = . (30) 

This means that the (forward-) difference scheme is symplectic. It can be proved that the scheme with respect to the 
(backward-) difference of q^"'\ 



A,aW=g« = l(oW-o('=-i)) 



(31) 
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is also symplectic as well. 

4 On the Symplectic Schemes 

We now show some well-known symplectic schemes in the Lagrangian formalism for the Hamiltonian systems. 

4.1 The Euler Mid-point Scheme for Separable Hamiltonian Systems 

The well-known Euler mid-point scheme for the separable Hamiltonian systems is as follows: 

h / -I- ni 



, = -hV,{^ ^) , (32) 

from which it follows that 

(n-l-l) _|_ ^(n) _|_ ~ 



fe'-^->-/-) + (/--p"-') = -/.F.( '' ) ' ), (33) 



Now it is easy to get the Euler mid-point scheme in the Lagrangian formalism, 



^ 2 )\ ■ (35) 

In order to show that it is symplectic, we first introduce the DEL one-form as follows: 

4") = dg-("){5("+^) - 2,(") + + \ [V,{^ + )] } . (36) 

Then the DEL condition gives rise to 

dg("+i) A dg^(") + dg("-i) A dg^(") = -:^yJ;+i/2) a dg^(") - A dg^(") , (37) 

where 

.("+1) _L ~ 



yin+1/2) ^y(rZ2±^y 



That is 

Now it is easy to prove that 

Ap^n+i) A dQ("+i) = (1 + ^l/(^"+i/2)^ d(?^("+i) A dg(") . (39) 
Therefore, the Euler mid-point scheme is symplectic. 

4.2 The Euler Mid-point Scheme for Generic Hamiltonian Systems 

For the general Hamiltonian H, the similar preserved symplectic form can also be given. Let us start with 

y p(ri+l) {n) „{n+l) _^ {n) 
qin+1) _ q(n) ^ ^^^^P 9 

/ . ^ /r,(»+i) _|_ r,(") n^n+l) _i_ Jn) . 

Introduce a pair of DEL one-forms 

4")(g) = d,^("){g("+i) - - hH,[^ ^ ^)} , 

^ N r /nC^+l) -I- n(") «("■+!) 4- \ ^ 

4"^ (P) = dp^'") - + hH, , ^-^^) } . (41) 

The DEL conditions for the pair of the DEL one-forms now read 
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From these conditions it follows that 
1 



dq 



(n+1) 



{dq^-'Y A [\hH^r'^'^ + 2(1 + lhH^;-y'^){hH(;-'/^Y'{^ - lhH(r,-'^'))] dg^ . (43) 



This shows that the following two- form in {dq^''^) is preserved, 



n) 



(dg("-^))^ A [^hH^r'^'^ + 2(1 + ^hH(;-y'^){hH(;-y'Y'(l - ^hH(r,-'/'^)] dq( 
It can be shown that it is nothing but the preserved symplectic structure, 

2(dp("+i))^ A d5("+i) = -{dq^^Y A [i/i£r(^+i/2) 

+ 2(1 + i/.i?(^+V2))(/,^(n+i/2))-i(i _ l/,if(«+V2))] dg("+i) , 
In terms of z"^ = {q^,p^), the mid-point scheme can be expressed as 

The DEL one-form for the scheme at tk now becomes 



(44) 



))}• 



(45) 



(46) 



(47) 



It is now straightforward to show that the symplectic structure-preserving law 

At(d0('=)^ A Jd^(^)) = 



holds if and only if the DEL form is closed. 

4.3 The High Order Symplectic Schemes 

Similarly, it can be checked that the high order symplectic schemes preserve also some two- forms in dq^''^ which 
are in fact the symplectic structures. Let us consider two examples for this point. 

The first one is proposed by Feng et al. in terms of generating function. 1^1 The scheme is as follows: 

^(n+i) = ^(") + /ij-iv,i?(i(^("+i) +^("))) - ^J-iV,((V,ir)^Jif,,JV,if)(i(z("+^) +^("))) . (48) 



In this case Ti. can be introduced as 

n = H- ^{V,HfjH,,JW,H . 
Then the fourth-order symplectic scheme can be rewritten as 

^(n+l) ^ ^(n) ^ /iJ-lV^H(i(^("+l) + ^("))) . 

Introducing an associated DEL form. 



.(fc) 



)^{ JAtzW - V,HW (i(^('=+^) + z^'^^)) } 



(49) 



(50) 



(51) 



it is easy to see that Ed\ and Ed2 differ an exact form. 



(52) 
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method is nothing but the mid-point scheme. Second, the stage two and order-four RK method is as foUows: 



Yi = + h 

Y2 = + h 

It can be expressed in terms of Hamiltonian H as 



=p(»)-i[H,{P„g,) + -ff,(P2,02)l, 



where 



Qi 


= -t 


-h 




Q2) 


Pi 


= p(")- 


h 




Q2) 




= H- 








Q2 


-h 




Q2) 


P2 


= p(")- 


h 




Q2) 



Introducing a pair of the DEL one-forms 

E^^\q) : = dp("){g("+i) - - ^{Hp{PuQi) + H^{P2,Q2))} 

E^S\p) : = d9("){p("+i) + ^{H,{P,,Q,) + H,{P2,Q2))} 

then the DEL conditions, i.e., their closed condition 

gives rise to the symplectic preserving property 

dp("+i) A dg("+i) = dp(") A dg(") . 



(53) 



(54) 



(55) 



(56) 



(57) 



(58) 



It can also be shown that 0;'"+^^ = dp^"'~^^^ A dg'^"+^^ may be expressed as dg'("+^^ A dg^"^ with some coefficients. 

5 The Trapezoidal Integrator 

It is well known that this scheme is good enough in comparison with other well-known symplectic schemes. But for 

the long time, it is not clear why it is so satisfactory. 

Wo will show that this scheme is symplectic, but the preserved symplectic structure is not simply ui = dp"'" A dq. 
Of course, the preserved symplectic structure should be canonically transformed to the one in the simple form with 
different canonical coordinates and momenta in principle. 

The scheme is given by 

q{n+l)_q{n) ^ ^[7J^"+l)(p(»+l),g(»+l)) + ijW(p(»),g(n))], 



(59) 



5.1 For Separable Hamiltonian Systems 
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the scheme reads 

qin+l) _ q{n) ^ + , p{n+l) _ p{n) ^ _ ^ (^(n+1) ^ ^ Vg^"\q^''^)] . (60) 

As what have been done before, let us introduce a pair of the EL one-forms, 

E^^\q) : = dp("){g("+i) - - , 

E^ip) : = dg("){p("+i) + - l^,^"^^^"^]} • (61) 

Then by some straightforward but more or less tedious calculations, it follows from the DEL condition, i.e. their closed 
property, that 

^(dg("))^(l + ^y^")) A (dg("+i) + dp("-^)) = -^(dg("))^(l + ^Vj^^) A (V^Jj^+i) dg("+i) + V}r'^dq<^-'M62) 
We get 

(dgW)^(l + ^fW) a (l + l^V^in+Dy^in+i) ^ _(d,W)^(i + ^^W) A (l + ^y^C^-D) d^^-^) . (63) 
This means that the following symplectic structure is preserved 

(d,("+i))^ (i + A (i + dg(") . (64) 

That is 

(dp("+i))^ A (l + ^4"+^)) d<7("+i) = (dp("))^ A (l + ^4")) dg(") . (65) 

It is straightforward to show that this two-form is closed and non-degenerate so that it is the preserved symplectic 
structure for this scheme. 

Using the following relation 

(dp("+i)) ^ (l + dg(«+i) - (l + dg(") , (66) 

one will get this two- form to be the same as Eq. (64). 

For the general separable Hamiltonian H = T{p) + V{q) we can get the preserved symplectic structure for the 
scheme as follows: 

= (dp("+i))^ A (i + :^rj;+i)v;(;+i)) dq^"+^^ . (67) 

It is also closed and nondegenerate. 

5.2 The Trapezoidal Scheme for General Hamiltonian Systems 

For the general system with non-separable Hamiltonian, the trapezoidal scheme gives 

q{n+l) _ q{n) ^ ^ jj{n)^^ _ ^(n) ^ _^(if + //(")) . (68) 

Similar to the separable Hamiltonian case, let us introduce a pair of DEL one-forms, 

E'^\q) : = dp("){g("+i) - - + <))} , 

E'^^^ip) : = d?W{p("+i) + ^(J?("+i) + Hf))} . (69) 

Similarly, by some straightforward but tedious calculations, the DEL condition for the pair of DEL one-forms gives 
rise to the following symplectic two-form and its preserving property. 



where 



= (dp("))^(l + ^H^;^H(^^ - '^H(,^^H(f - H^fH^j;^) A dg^ - ^(dg("))^ffWif(^) A dg^ . (71) 
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If we introduce new variables 

p(")=p(«)-^i/(«), 

it follows that 

wi"^ = dp(") A dg(") . (73) 

This is another expression for the preserved symplectic 
structure in the trapezoidal scheme. 

6 Some Remarks 

i) In order to show whether a scheme for a given 
Hamiltonian system is symplectic preserving, the 
first issue in our approach to be considered is to 
release the scheme from the solution space to the 
function space. Otherwise, it is difficult to make 
precise sense for the differential calculation in the 
solution space. One of the roles played by the 
EL cohomological concepts is just to release the 
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